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The generalized harmonic equations of general relativity are written in 3+1 form. The result is a 
system of partial differential equations with first order time and second order space derivatives for 
the spatial metric, extrinsic curvature, lapse function and shift vector, plus fields that represent the 
time derivatives of the lapse and shift. This allows for a direct comparison between the generalized 
harmonic and the Arnowitt-Deser-Misner formulations. The 3+1 generalized harmonic equations 
are also written in terms of conformal variables and compared to the Baumgarte-Shapiro-Shibata- 
Nakamura equations with moving puncture gauge conditions. 
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I. INTRODUCTION 

The generalized harmonic equations [J-Q are a sym- 
metric hyperbolic formulation of general relativity. They 
were originally written as a second order system of partial 
differential equations for the spacetime metric '"'g^ii/. By 
adding extra variables to represent derivatives of '""'g^iy, 
the generalized harmonic equations can be written as a 
fully first order system 0, Q , or as a system with first 
order time and second order space derivatives @. Typ- 
ically the fundamental variables are the components of 
the spacetime metric and its derivatives. 

In this paper we carry out a 3+1 splitting of the gen- 
eralized harmonic (GH) equations. In this way the GH 
system is written in terms of traditional 3+1 variables 
with first-order time and second-order space derivatives. 
The 3+1 variables include the spatial metric gij, extrinsic 
curvature Kij, lapse function a and shift vector /3*. The 
extrinsic curvature is directly related to the time deriva- 
tive of the spatial metric; likewise, we introduce fields tt 
and that are directly related to the time derivatives 
of a and . The result of this analysis is a concise and 
elegant expression of Einstein's theory. 

Currently there are two formulations of the Einstein 
equations in widespread use in the numerical relativity 
community. One is the generalized harmonic system, 
the other is the Baumgarte-Shapiro-Shibata-Nakamura 
(BSSN) system 0, i] along with moving puncture gauge 
conditions The BSSN equations are direct descen- 

dants of the Arnowitt-Deser-Misner (ADM) equations, 
which are obtained from a 3+ 1 sp litting of the Einstein 
equations [HI- (See also Refs. [ll[ll.) ADM and BSSN 
are typically written as systems with first-order time and 
second-order space derivatives. The fundamental vari- 
ables for ADM are the 3+1 variables gij, Kij, a and 
/3V BSSN is obtained from a change of variables, defined 
by conformal splitting, and the introduction of new in- 
dependent variables, namely, the conformal connection 
functions. Most often the BSSN system is supplemented 
with the moving puncture gauge conditions which take 
the form of evolution equations for the lapse function and 
shift vector. 

In earlier work, Friedrich and Rendall [l^l (see also 
Ref. [l^) wrote the generalized harmonic equations in 
terms of 3+1 variables gij, a and Their motiva- 



tion was not to compare GH to ADM or BSSN. Con- 
sequently, the relationship between the GH and ADM 
or BSSN systems has remained obscure. In Sec. HI the 
precise relationship between the GH equations and the 
ADM equations is presented. The relationship between 
the GH equations and the BSSN equations is displayed 
explicitly in Sec. V. Note that in Refs. [H, [T^ the Z4 
formulation [l^ of general relativity is written in 3+1 
form with a conformal splitting, and used to compare Z4 
to BSSN. 

In See. 2 we review the generalized harmonic formu- 
lation of general relativity and discuss its interpretation 
as an initial value problem. In Sec. 3 we write the GH 
equations in 3+1 form, compare the results to ADM, and 
show that the system is symmetric hyperbolic. Technical 
details are contained in Appendix A. In Sec. 4 the 3+1 
GH equations are written in terms of conformal variables. 
The GH equations are compared to BSSN and the mov- 
ing puncture gauge in Sec. 5. In Appendix B we show 
that the GH system with moving puncture gauge condi- 
tions has the same level of hyperbolicity as BSSN with 
the moving puncture gauge. A brief summary is provided 
in Sec. 6. 



II. GENERALIZED HARMONIC EQUATIONS 

Let '""'(^pi/ denote the physical spacetime metric with 
Christoffel symbols '^'T^o-p, covariant derivative V^, and 
Ricci curvature '-'^^R^^u- Let '"'f ''jrp denote a background 
connection; this connection might be built from a back- 
ground metric '"'^^y. As discussed in Ref. [3, the back- 
ground connection is needed for general covariance. For 
practical applications it would be natural to choose F'^g-p 
to be the flat connection. If the coordinates are inter- 
preted as Cartesian, then the components F^o-p are zero. 

Now introduce a spacetime vector field H^^, the "gauge 
source vector" , and define 



C^'F^<,p - '^'f ^p),9 



crp 



(1) 



Note that the physical and background connections only 
appear as the difference '■''F'^g-p 
forms as a tensor. 



'"T^o^p, which trans- 
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The generalized harmonic equations are 

+8^G[T^,-'*'5^,T;/2] , (2a) 
= . (2b) 

The term proportional to Newton's constant G repre- 
sents the matter content, where T^i^ is the matter stress- 
energy-momentum tensor. The matter equations of mo- 
tion imply the conservation laws V ^T^^ = 0. The term 
proportional to the constant k enforces constraint damp- 
ing; it depends on a timelikc, future-pointing unit vector 
field . Below we assume that this vector field is the 
unit normal to a set of spacelike hypersurfaces t = const. 

The GH equations ^ are equivalent to the Einstein 
equations = SttG [T^^ ~ <*'g^^r^/2]. This follows 

trivially by inserting Eq. ^p) into Eq. ((2^). What makes 
the GH equations interesting, and useful, is their inter- 
pretation as an initial value problem. 

Let us define M^, = -('">G^^ -87rT^i.)n'' where '^'G^^ 
is the Einstein tensor and is the unit normal to the 
t = const slices. Note that H = —2M^n^ and Mi are the 
Hamiltonian and momentum constraints, respectively. 

The initial value interpretation of the GH equations 
relies on two key results. The first is obtained by con- 
tracting Eq. (I2K) with the unit normal ti^. This yields 
an equation of the form [l^ 

^tC^' = {terms - X, C, 9,C} (3) 

where the terms on the right-hand side are proportional 
to the constraints and C^, and spatial derivatives 
of . This equation can be rearranged to show that 
M.^^ is proportional to C^, the time derivative of C^, and 
spatial derivatives of C^. It follows that for any solution 
of the GH equations ([2]), the Hamiltonian and momentum 
constraints A^'' = must hold. 

The second key result is obtained from the covariant 
derivative of Eq. . After applying the Ricci identity 
and using the result ([3|) , we find [a, [l^ 

dtM^" = {terms - M, d^M, C, diC, d^dfi] . (4) 

Together, Eqs. ^ and @ show that as long as the con- 
straints = and A^^ = hold at the initial time, 
then they will continue to hold for all time. 

From the preceding analysis we see that a solution of 
Einstein's equations can be found by choosing initial data 
that satisfy both sets of constraints, C = and M.^ ~ 0, 
at the initial time, and then evolving this data into the 
future using the GH equation ([2^). 

The generalized harmonic formulation of general rela- 
tivity is important because the GH equations are sym- 
metric hyperbolic, provided the gauge source vector i?^ 
is specified directly as a function of the spacetime coor- 
dinates and metric ''''g^i/. In particular, the second 
derivative terms in Eq. ((2^) combine to form a wave op- 
erator '•'^^ g'^'^dcrdp acting on the spacetime metric '"'g^iy. 
If, on the other hand, the i/^'s are specified directly and 



depend on d^'^'^^ g^i,^, then Eq. ((2^) will include terms that 
interfere with the nice wave operator. In general the sys- 
tem will no longer be symmetric hyperbolic. 

In much of the early numerical work with the GH equa- 
tions, iJ^ was not specified directly. Rather, was ele- 
vated to the status of a dynamical variable by introducing 
"driver" equations. With a driver equation, V^V^i/" or 
dtH'^ is set equal to some function of '^'g^^i/ and and 
their derivatives [1, [l^l • In this case the analysis of hy- 
perbolicity is more complicated. 

Recent work by Szilagyi. Lindblom and Scheel [2l| has 
shown the practical benefits of the "damped wave gauge" . 
For this gauge condition is specified directly as a 
function of the spacetime metric. Throughout this paper 
I will assume that the gauge source vector is specified 
directly. If it depends only on the coordinates x'^ and 
metric '"'g^i/, then the system is symmetric hyperbolic. 
In Sec. 5 we consider the GH equations with the moving 
puncture gauge. In this case i/^ depends on derivatives 
of the metric, and the system is not symmetric hyper- 
bolic. In Appendix B we show that this system has the 
same level of hyperbolicity as BSSN with the moving 
puncture gauge. 



III. GH EQUATIONS IN 3+1 FORM 

Let us begin by reviewing the 3+1 decomposition of the 
Einstein equations pH - fTsj . The analysis yields evolution 
equations 

d±gij = -2aKij , (5a) 
d^K.j = a [R,j - 2K,kK'y + KK,j\ - D.D^a 

-8^Ga[sy -g,,(s-p)/2] (5b) 

for the spatial metric gij and extrinsic curvature Kij. 
Here, Di and Rij denote the covariant derivative and 
Ricci tensor built from the spatial metric. The lapse 
function is a and the shift vector is The time deriva- 
tive operator is defined hy d± = dt — Cp^ where £^3 is 
the Lie derivative along the shift. The matter variables 
are the energy density p = n^n^T^j^y^ momentum density 
ji = —n^^T^i, and spatial stress Sy = Ty. The 3+1 split- 
ting of the matter conservation equations V^T'"^ ~ 
gives [111 

d^p = as^^ Kij + apK — aDif — 2fDia , (6a) 
d^ji — otKji — sijD^ a — pDiU — aD^ Sij . (6b) 

The spatial metric and extrinsic curvature must also sat- 
isfy the Hamiltonian and momentum constraints, 

H = K^ -K,jK'^ +R~16ttGp = , (7a) 
M., = DjKl - D,K~87rGj, = . (7b) 

If the constraints hold at the initial time, then the evolu- 
tion equations ^ and © insure that they will continue 
to hold at future times. 
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In the numerical relativity community the results ([5]) 
are referred to as the Arnowitt-Deser-Misner (ADM) 
equations Here we use the common convention of 

writing these equations in the form used by Smarr and 
York [il,!!!. 

The mathematical details of the 3+1 splitting of the 
generalized harmonic equations ^ are presented in Ap- 
pendix A. The result is the following system of evolution 
equations, 



-2a A' 



d±K.,j = a [R^j - 2K^kKj - nluj] - D^Dj 
—aD(^iCj) —KagijC± /2 
-87rGa [sij - gij{s - p)/2] , 



(8a) 



(8b) 
(8c) 
(8d) 



d±a 

dtP' = fi^Djfi' + a^p' - aD'a + a^H' 
d±T: = -aK.jK'^ + D^D'a + CD,a 

-KaC±/2-ATrGa{p + s) , (8e) 
d±p' = g^^DuDiP' + aD'TT - irD'a - 2K^Wja 



-2aK^''Arjk+KaC'-167TGaf 



and constraints, 



C = -p' + AY'jug''' 



R-16ttGp , 



Mi = DjKl - D^K^8^^Gj^ 



(8f) 



(9a) 
(9b) 
(9c) 
(9d) 



The dependent variables include the spatial metric g^ , 
extrinsic curvature Kij, lapse function a and shift vec- 
tor /3*. Wc have also introduced the variables tt and p*. 
Equation ^) shows that tt is related to the time deriva- 
tive of a. Likewise, from Eq. ([5Ji) we see that is related 
to the time derivative of /3*. Note that the gauge source 
vector H'^ appears in these equations as a spatial scalar 
H± and a spatial vector The source H± appears 

in the evolution equation ([5J;) for the lapse a, while the 
source appears in the evolution equation ((Sji) for the 
shift 13\ 

In deriving the 3+1 GH equations ([5]), (jH]) we have 
assumed that the only non vanishing components of 
the background connection '""F^o-p are the spatial com- 
ponents ''''f'jfc. This is equivalent to building the 
background connection from a background metric ^'^^g^u 
which, under a 3+1 splitting, has unit lapse, vanishing 
shift, and a time-independent spatial metric. In this case 
the only remaining background structure is the spatial 
connection whose components are f = '"''f '^fc. We also 
assume that the background spatial connection is flat, 
and in Eqs. (|5]), ^ use the notation 



jk 



(10) 



Finally, we let Di denote the covariant derivative built 
from the background connection. 



Comparing the 3+1 GH equations ([8^) and ^p) with 
the ADM equations ([5]), we find 



id±9i])GH - id±gij)ADM = , 

ADM 



(11a) 



'Kag,fi^l2 . (lib) 



As expected, the difference is proportional to the con- 
straints dH]). 

The constraint evolution system for the 3+1 general- 
ized harmonic equations is 



9_lC_l = -aA'Ci + a-H + CD,a - aDiC 
--2KaC± , 

d±C,, = C±D,a - aD,C± - 2aMi - 2aKijC' 



(12a) 
(12b) 



—KaCi , 

d^n = -2a7:n + 2aRCi_ - iM^D'a 

-2aD'M^ + 2a{ie^ - Kg'^)DiCj 
~2KaC±-32TTGapC± , (12c) 



d±M., = -nD,a + (KSj - Kl)Dj{aCA 



1 



-aDM 



-anM, + D^aD[,Cj] + Di[aDp) 
-KDi{aC±)—8TrGajiC± . 



(12d) 



These results are found from the evolution equations ([5]) 
and ([6|) applied to the definitions ([9]). 

The GH equations are symmetric hyperbolic. This can 
be shown by considering the second-order system ([2^), 
or the fully first order system of Refs. 0, [E]- Gundlach 
and Martfn-Garcfa |23| have given a definition of sym- 
metric hyperbolicity that applies to quasilinear systems 
of partial differential equations with first-order time and 
second-order space derivatives. We can apply their defi- 
nition to the 3+1 GH equations ([5]). 

To begin, we assume that the matter fields are not 
derivatively coupled to gravity; that is, the matter La- 
grangian does not contain derivatives of the metric. Then 
the matter variables p, j', and Sij do not contain deriva- 
tives of the gravitational variables gij, Kij, a, tt, /?% and 
p'. We also assume, as discussed in Sec. 2, that the gauge 
sources and H± arc directly specified in terms of the 
spacetime coordinates and the metric variables gij, a, 
and /3*, not on their derivatives. 

The analysis can be described as follows. In effect, we 
assign weight to the metric variables and weight 1 to 
the "velocities" Kij, tt, and p'. One unit of weight is 
added for each derivative. We introduce the weight 1 
variables 



gmij — (^mgij i 

ai = diU , 



(13a) 
(13b) 
(13c) 



defined as derivatives of the weight variables, and 
compute their equations of motion by differentiating 
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Eqs. (H^), (Ht) and ([8ji). Note that diUj, dig^ke, and 
di{l3jkg''^) are symmetric in i and j. The principal parts 
of the GH equations ^ arc constructed from the highest 
weight terms in the equations of motion for the weight 1 
variables; these are 



dtiT 



29„ 



2admKij 



2 

a 



adiaj 



dtPi = adiTT + g^^djjSut 



(14a) 

(14b) 

(14c) 
(14d) 
(14c) 
(14f) 



where = denotes equality apart from lower weight terms. 
Here we have defined the operator dt = dt — P^dk- 
We now define the quadratic form 



kl 



~M 



-ig^'^gmij 9nM 



-^g Pki Pij 

+ {pi - ai/a){pj - aj/a) 



(15) 



where the tensors M*-'*^^, A'/*^ and M (not related to one 
another) are positive definite. Direct calculation using 
Eqs. p4)) shows that the time derivative of e has a prin- 
cipal part that can be written as the gradient of a vector 
0'; that is, dt£ = di<j)^ . This shows that e is a quadratic, 



positive-definite energy density with flux <jf . It follows 
from the theorems of Gundlach and Martm-Garcia [l^] 
that the system ([8]) is symmetric hyperbolic. 



IV. GH EQUATIONS IN CONFORMAL 
VARIABLES 

In 3+1 form the GH variables are , Kij, a, tt, p'. 
Introduce a time-independent spatial density of weight 
2, denoted 7. As this notation suggests, 7 can be cho- 
sen as the determinant of a background metric 7^^ , and 
this same background metric can be used to define the 
background connection V^ij- Now consider the confor- 
mal variables 7^^, Aij, Lp, K, A*, a, tt, and /?' defined by 



A, 



ii/g) 



1/3 



K, 



;9ijK 



(16a) 
(16b) 

(16c) 



" 12 ' 
A' = {g/l)'^'p' + l{g/j)-'^'D\g/j) . (16d) 

Note that the determinant of jij is 7 and the trace of Aij 
vanishes. The inverse relations are 



9ij 



s^^ A 



(17a) 
(17b) 
(17c) 



Indices on the new variables Aij and A* arc raised and 
lowered with the conformal metric 7^ j . 

In terms of the conformal variables, the GH equations 
(El) arc 



d 



d^A, 



6 6 
y 3 



D^a + 2D'^D,a + a (n - KCi_ - D.C - 6C*9,(^) 



■ 3aKCi/2 -h 47rGa(/9 + s) , 



n TF 



aTZij — 2aDiDjip + AaDi(pDj(p — DiDjU + ADi^aDj\(p — SnGas 



2 ~ - , ~ -I, 

--A.jDkP^ - 2aA,kA'] + aKA,^ - aC±A,j + , 
o 



4C(,^,)¥>-CfcAf(, 



TF 



5j.A' = 7"5fe5«/3' + l^'Dkfi'' + ^D\Dk/3^) - 2A'''dka + 2ai'=^Af + UaA^'^dk^p 

-^aD'K + aD'C± + ^ae^'^KC + nae^^C - IQnGae^-^f , 
d±a = a^TT — a^H± , 



(18a) 
(18b) 

(18c) 
(18d) 

(18e) 
(18f) 
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where 



3Wi 



d±7T ^ -aAiiA'^ - - 
u 3 



aK'^ + e 



A' - 2L)V - D'a/a , (18g) 
(^D^a + 2£)VAa) + CD,a - KaC±/2 - AnGaip + s) , (18h) 



,j _ --r'DkD,%j + 7*'^[Ar"fe,Ar(,,)„, + 2Ar'"fe(,Ar,)„, + Ar^.^Ar^j,] + ^u(,d^)K^ . (19) 



r 



Wc have also defined 



rl Ti? 



where FJ^. are the Christoffel symbols built from the con- 
formal metric. 



where R is the Ricci scalar built from the conformal met- 
ric 7ij. One must remember that the indices on Ci and 
(20) J\Ai are raised and lowered with the physical metric gij. 
Thus, for example, C' = e-'''^(-A* + l^P^kl^''). 

The system IjlSp is, of course, symmetric hyperbolic 
as long as H±_ and Hi do not depend on the weight 1 



In terms of conformal variables, the constraints (|9]) are variables {K , Aij, A*, tt) or derivatives of the weight 

variables {"fij, ip, a, We can confirm this by defining 



C±=TT + K , 

n = - A,jA'^ - WnGp 

o 

2 ~ 



(21a) 
(21b) 

(21c) 



ai = dia , 
<Pi = diip , 



(22a) 
(22b) 
(22c) 
(22d) 



and computing the principal parts of the evolutions equa- 



Mi = DjA\ - -DiK + QAlDjip - SnGji , (21d) tions for the weight 1 variables: 



dt^Pn 

u 

dtK ae-^ff^ 



^f^dmPij ~ \admK 
6 6 



9iAj — 89i(y5j dia 



dtA, 



,1 



-ae '^'^-i^^dkltij + ae 



a 



^^Aj - 2di(pj diUj 

a 



TF 



4A, = 7"afe/3« + i7"5,/3fe£ - ^ad^k + ad^ir , 



2 -4m 

a e ^ 



d,Kj - 2d, 



ipj diUj 

a 



One can show by explicit calculation that the positive definite energy density 



(23a) 
(23b) 
(23c) 

(23d) 

(23e) 
(23f) 
(23g) 
(23h) 



+M 



^7""(7mij + 4%-(/7™)(7„fef + 47M(/9„) + e'^'^{A,j + ^ijK/3 - /3(y)/a)(ifc£ + %iK/3 - /3(fc£)/a) 



^e^'^7"/3fc,/3f, + (A, - 2ip, - a,/a)(A, - 2(^, - a, /a) 



in)' 



7 "iaj 



(24) 



satisfies the conservation equation dtS = dii 



V. COMPARISON WITH BSSN AND THE 
MOVING PUNCTURE GAUGE 



The BSSN equations in covariant form are [1^ 



6 



d±K 



6 6 



1 



D^a + 2D'ipD,a 



+ 47rGa(/9 + s) 



-iip 



alZij — 2aDiDj(p + AaDiipDjip — DiDja + 4:D(^iaDj^ip — SirGaSi 



TF 



^-A.jDkP" - 2aA,kA^ + aKA.j , 



+2aA'"'\V'M + l2aA'^dkip - -aD' K - IGnGae^"^ f 



(25a) 
(25b) 
(25c) 

(25d) 
(25e) 



The variables A* are the "conformal connection func- 
tions" . If the background is flat and the coordinates are 
interpreted as Cartesian, then the background connec- 
tion vanishes, T'^jk = 0. (We also have Di — di.) In this 
case it is common to use the notation F* = jkl^^ for 
these variables rather than A' . Also observe that the first 
term on the right-hand side of Eq. (HS^), —j^'^CjDkp'^ ~ 
AWj(3'' - j'^^AT^kiDjP'', and the Lie derivative term on 
the left-hand side, -CpA' = -pWjA' + AWjl3\ com- 
bine to insure that only derivatives of A*, and not A' 
itself, appear in Eq. (PSb ). This rule is discussed in 
Ref. [10[ and is followed by most numerical relativity 
groups who use the BSSN system. 

The BSSN equations are usually accompanied by the 
moving puncture gauge conditions. 



dta = P'd.a - 2aK 



(26a) 
(26b) 



I 

where 77 is a parameter, independent of the field vari- 
ables. Eqs. p6k) and ([26}d) are the 1-l-log slicing 
and the gamma-driver shift conditions, respectively. The 
gamma-driver shift is often written as a system of first- 
order equations for the shift vector /3* and an auxiliary 
field As shown in Ref. [2^, these equations can 

be integrated to yield the single equation (|26b ') for /3*. 



By explicitly comparing the GH equations in conformal 
variables, Eqs. (fT8|) . with the BSSN equations ([25]), we 
find 



{daij)GH - 


- (pi'^ij) BSSN 


= , 




- {dt(p)BSSN 


= , 




- {dtK)BSSN 


^a{H- KCi_ 


{dtAij)cH ~ 


{dtAij)BSSN 


= —aCj_Aij + 


{dt~A')oH - 




= f%Dkf3' - 



TF 



(27a) 
(27b) 
(27c) 

(27d) 
(27e) 



As expected, the differences between GH and BSSN are 
proportional to the constraints. Note that the terms 
proportional to C± simply exchange tt for —K; like- 
wise, the terms proportional to Ci simply exchange A' 
for ATijkj-''^- Also observe that only a few of the terms 
on the right-hand sides of Eqs. (P?]) contribute to the 



principal parts of the equations. In particular, we have 

idtK)aH - idtK)BssN = a{n - AC') (28a) 

{dtA')aH - idtA')BssN = a&Ci_ . (28b) 

The principal parts of the GH and BSSN equations for 
7y, (p, and Aij coincide. 

The results (j27p provide a simple prescription for con- 
verting a BSSN code into a GH code. First, add the 
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terms on the right-hand sides of Eqs. (|27p to the BSSN 
equations of motion. Next, add the equation of motion 
(|18h ) for TT. Finally, modify the evolution equations for 
a and so that they take the form of Eqs. (HHf) and 

With an appropriate choice of the gauge sources H± 
and Hi, we can adopt moving puncture gauge conditions 
within the generalized harmonic formalismjj In terms of 
conformal variables, we need 



H^=TT + 2K/a 



4a2 



A' 



V 



13' 



(29a) 



(29b) 



so that the GH equations (HHf) and ([T5k ) coincide with 
the moving puncture equations In terms of the 

original 3+1 variables, we have 



Hj_ = TT + 2K/a 
3 

4q,2 



{9/iy^'p' + li9/ir'^'D\g/j) 



-p' + D^a/a 



and the moving puncture gauge conditions read 
dta = P'd^a - 2aK , 

3 

+ 4 



{9/iy^'p' + li9/l)-'^'D\g/^) 



(30a) 

(30b) 

(31a) 
.(31b) 



With the moving puncture gauge, the H's depend on 
weight 1 variables and derivatives of weight vari- 
ables. This spoils the symmetric hyperbolicity of the 
system. In Appendix B we analyze the GH equations 
with moving puncture gauge conditions and show that 
they are strongly hyperbolic as long as the condition 
2a ^ {9/lY^^ is met. Note that one can use the con- 
straint = to exchange for Ar^jfegJ*^ in Eq. (pib). 
This does not affect the hyperbolicity of the system. 



VI. SUMMARY 

The generalized harmonic equations have been written 
in 3-1-1 form using as independent variables the spatial 
metric gij , extrinsic curvature K^j , lapse function a and 
shift vector as well as fields tt and related to the 
time derivatives of a and The resulting set of evolu- 
tion equations ([8]) and constraints ^ are a concise and 
elegant formulation of general relativity. The GH evolu- 
tion system is symmetric hyperbolic with the conserved, 
positive definite energy density displayed in Eq. ([15]). 



^ Moving puncture gauge conditions for the Z4 formulation of gen- 
eral relativity have been discussed in Refs. [Tsl. llTl. [27|. 



The 3-1-1 GH equations are written in terms of confor- 
mal variables in Eqs. (fT8|) and ([21]) . This allows for a di- 
rect comparison with the BSSN formulation of Einstein's 
theory, and provides a simple prescription for convert- 
ing a BSSN code into a GH code. The moving puncture 
gauge conditions cannot be used with the GH equations 
without spoiling symmetric hyperbolicity. Nevertheless, 
the GH system with moving puncture gauge has the same 
level of hyperbolicity as the BSSN system with moving 
puncture gauge. 
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Appendix A: 3+1 splitting 

In this appendix we derive the equations of motion ([5]) 
and the constraints ([9]) by carrying out a 3+1 splitting 
of the spacetime generalized harmonic equations ([5]). 

Let = — Qf5^ denote the covariant normal to the 
spacelike hypersurfaces t = const; the contravariant form 
is = ((5^ — /3*(5f)/a. Also introduce the operator 
= 5^ that projects spacetime covectors into spa- 
tial covectors. The covariant form of this operator is 

0. 



<5; + /3M, 



,,, it satisfies X'lX^^ 



5) and X],n^' 

The spacetime metric is written in terms of the normal. 



spatial projection operator, and spatial metric as 



(Al) 



Spacetime indices p,, v, etc. are always raised and lowered 
with the spacetime metric '"'g^i/ and its inverse '"'g'"^, 
while spatial indices i, j, etc. are always raised and low- 
ered with the spatial metric g^j and its inverse g*-' . 

The spacetime Christoffel symbols can be written in 
terms of 3+1 quantities as 



(A2a) 

(A2b) 
(A2c) 

(A2d) 
(A2e) 

x;<^'r^.p7i-x; = -a- + id,n/a , (A2f) 

where d± = dt—Cp is the time derivative operator used in 
the main text. The results (|A2[) are obtained by comput- 
ing the normal and tangential projections of derivatives 
of the spacetime metric, 9(7'"' 5^1/, and using the defini- 
tion of the Christoffel symbols. Also note that we have 
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used the relation d±gij = —2aKij that defines the ex- 
trinsie curvature. This is the equation of motion (|5^) for 
the spatial metric. 

We will also need the splitting of the Ricci tensor, 

'■^^Rf.^nf'n'' = {d^K + D,D'a)/a 

-K,jK'^ , (A3a) 

-{d^K,.j)/a - {D,D,a)/a , (A3b) 
<^'i?^,Xf n'' = -DjKI + D,K , (A3c) 

and the curvature scalar: 

<^'i? = i? + K\jK'^ - 2{d^K)/a 

-2{D,D'a)/a . (A4) 

These results can be obtained from the definition of 
the Riemann tensor in terms of covariant derivatives, 
^^^Rf^^apVP = V^V^K - V^^V^K, or from the defini- 
tion of Riemann in terms of Christoffel symbols and the 
resuhs (|X2|) . 

The GH constraint C^. = iJ^ + ('^T^^p - '"'f 
must be split into a normal component, C± = C^?!^, and 
a tangential component, Ci = C^X^ . These calculations 
depend on the 3+1 splitting of the background connec- 



tion T'^ap. Let us assume that the background connec- 
tion is constructed from a background metric '^'g^i/. This 
metric can be split with respect to the t = const hy- 
persurfaces into the 3-1-1 quantities gij, a, and /3\ The 
results (|A2[) . applied to the background geometry, can 
be rearranged to give the components of the background 
connection: 



'^>f = {dta + P'd,a - P'B'K,j)/a , (A5a) 

<^'f = {d,a - P'k,,)/a , (A5b) 

<^'f *y = ~K^J/a , (A5c) 
'■'^T\t=a&a-2aP^Kjkg''' 

+dt$'+pWjP\ (A5d) 

-akjkt' + DjB' , (A5e) 

<^'f%- =f^, . (A5f) 



Here, the background extrinsic curvature is defined by 
{dt - Cp)gij = -2aK^j. 

The calculations for the normal and tangential compo- 
nents of the constraint yield 



C±=H± + K + X;d^a - '^g'^K,, -^{dt- — A^^A/3^\fC„ ^A/S^a^a , (A6a) 

OL aa aa aa 

a a a a 

+ 4^3,, A/3^" [dta - {2(3'' - p'')dka + a^g^^'Ru] - 4^5.^ A/?'^ [A/3^ A/3^ - 2o?g^'] Km ■ (A6b) 



Here we have defined iJ_L = i/^rip, and A/?' = /?' — /3'. 

Note that each term in Eq. (|A6b ,) is a spatial scalar, 
and each term in Eq. (jA6b ) is a spatial covector. In these 
equations we can absorb terms that depend on the phys- 
ical tensors 5^, a, the background tensors gij, a, /3^, 
and derivatives of these background tensors into H± and 
Hi. We cannot absorb terms that depend on derivatives 
of g.ij, a or /3* because this would change the hyperbol- 
icity of the GH system. Thus, we have the following 
results: 



I 

Let us define 

TT = ^d±a + H± , (A8a) 

p, = ^g,,{dt(3' ~l3''D,p^) 

+ -^^a - H, . (A8b) 
a 

When rearranged, these definitions become the equations 
of motion ([5J;) and ([Hli) for a and /3'. The constraints 
become 



1 



C± = TT + A' , (A9a) 
Cx = F± + A- + —d^a , (A7a) c, = -p, + Ar,,,.9^'= , (A9b) 

C, = H, + AFy-feg^*'- - -d,a 
a 

'-^dijidtP^ - P'^DkP^) . (A7b) Our next task is to split the terms V(^C,,). The normal 



which are Eqs. (jH^) and ([Hb) from the main text. 
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and tangential projections are 



-idi_C^-Cdia) , (AlOa) 
a 



1 

1 

+ 2 



(AlOb) 

a^Ci + i^TijC^' . (AlOc) 



The projections of the spacetime GH equation ((2^) are 
obtained from the Eqs. (|X9)) and (jA10|) above. The 

result for the normal-normal projection is 

-KaCi/2-47rGa(/9 + s) , (All) 

which is Eq. ([HJj) from the main text. The tangential- 
tangential projection yields 

-aD( fij-^ - Kacjt J Cx /2 



D.,Dja 



-87rGa [s^j - (s - /3)/2] 



(A12) 



which is Eq. (|Hb)- 

The normal-tangential projection of the spacetime GH 
equation leads to the result 

d^p' = g'^'^DkDeAlS' + aD'n - nD'a - 2K'Wja 

-A/3"i?„fe,,] . (A13) 

We now assume the background lapse is unity, a = 1, 
and the background shift vanishes, = 0. We also as- 
sume that the background spatial metric gij is flat and 
time independent. These assumptions imply that the 
background extrinsic curvature Kij and background Rie- 
mann tensor RmMj vanish. Then the normal-tangential 
projection becomes 

d±p' = g'^'^DkDiP' + aD'n - nD'a - 2K'Wja 

+2aK^''Arjk+KaC'-16nGaf , (A14) 

which is Eq. (|8f) from the main text. 

The analysis shows that the spacetime GH equations 
^ are equivalent to the evolution equations ([8]) plus the 
constraints Cj_ = and = 0. The constraint evolution 
system ([9]) shows that C± ~ and Ci = will hold for 
all time if and only if all of the constraint functions C± , 
Ci, Ji, and Mi vanish. It is sufficient to impose these 
constraints at the initial time; the evolution equations 
will insure that they continue to hold into the future. 



Appendix B: Hyperbolicity of the GH equations 
with moving puncture gauge 

In this section we analyze the hyperbolicity of the gen- 
eralized harmonic equations with the moving puncture 



gauge conditions ((3T|) . That is, we consider Eqs. ((8)) with 
the gauge sources H±^ and Hi given by Eqs. Sym- 
metric hyperbolicity is spoiled by the presence of tt, Ktj, 
and derivatives of gij and a in the 's. Nevertheless, 
the equations form a quasilinear system of partial differ- 
ential equations with first order time and second order 
space derivatives. We can ap ply the pseudo-differential 
reduction techniques of Refs. [2^-[30l to check for strong 
hyperbolicity. 

The principal parts of the equations are constructed 
from the highest weight terms. We identify the "coordi- 
nate variables" gij , a and as weight and the "veloc- 
ity variables" Kij, tt and as weight 1. Each derivative 
adds a unit of weight. The principal parts of the GH 
equations with moving puncture gauge conditions are 



dtgij = 2gk(^,dj)l3'^ - 2aK., 



I] , 



(Bla) 



dta 



a 



g dkdigij + ad(iPj) - did^a , (Bib) 



2 

-2a A' , 



P' + -^g^'g^'d.gk, 



Btn = g^^didja , 
dtp' = ag'^d,n + g^''d,dkP' 



(Blc) 

(Bid) 

(Ble) 
(Blf) 



where dt = dt — P'^di. 

Let Hi denote a covector normalized by the spatial met- 
ric: Hig^^rij = 1. The principal symbol for the system 
(|B1[) above is defined by 



(B2a) 



fiKij 



;9t] + 0!gk{znj)p'' - n^rija , (B2b) 



-2ag'^K, 



= ^(,9/7)^/^ 



P" + -^n'g''^gkt 



(B2c) 

(B2d) 

(B2e) 
(B2f) 



where jl = p — /3^ni. The proper speed (proper distance 
per unit proper time as measured by observers at rest in 
the t = const slices) of a characteristic mode is given 
by iP'^ — p)/a. (See, for example, the discussion in 

Ref. fai). 



Now introduce an ortlionormal triad consisting of rii 
and unit vectors e\, with A = 1,2. These vectors sat- 
isfy nie\ = and eAdv^B ~ ^ab- When we project 
equations (jB2p into the triad directions and e\, the 
principal symbol separates into blocks that have com- 
mon transformation properties under rotations about the 
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plane orthogonal to . The scalar block is 





= 2/3^ - 2aK^i_ , 


(B3a) 




= -2akABS^'' , 


(B3b) 






(B3c) 




^ ,-4 R 

= -^gABS^"" , 


(B3d) 


[lot 


= -2a{k^^ + kABS^'') , 


(B3e) 










+l{g/i)'^Hg±± + gABS^^) 


, (B3f) 


/iTT 


= a , 


(B3g) 




= a7r + ^-^ . 


(B3h) 



Here and below, the _L and upper case Latin indices 
are defined, for example, by g±± = gijU^n^ and cjab = 
gije\e!'Q. The vector block is 



IJ-gi-A = Pa - 2aK±A , 
a a ^ 

= -^g±A + ^Pa , 

f^^A - l{g/i)'^'PA , 



PPA = Pa 
The tensor block is 



uk*f ---f/f 

P-^AB — 2^^B 



AB 1 



(B4a) 
(B4b) 

(B4c) 
(B4d) 



(B5a) 
(B5b) 



where the superscript tf indicates that the tensor is 
trace-free in the two-dimensional surface orthogonal to 

Hi. 

A quasilinear system is strongly hyperbolic if its princi- 
pal symbol possesses a complete set of eigenvectors with 
real eigenvalues /i. The tensor block (|B5|) meets these 



criteria with eigenvalues fi = [3-^ ± a. These eigen- 
values correspond to proper speeds of ±1. The vec- 
tor block also meets the criteria for strong hypcrbolic- 
ity with /i = ± a and /i = ± \/3(g/7)i/3/2. 
The proper speeds for the vector modes are ±1 and 
±V3U/W^/{2a). 

The eigenvalues for the scalar block arc fi = ± a 
(with multiplicity two), fj, = ± ^/2a, and /i = (5^ ± 
{g/lY^^- These correspond to proper speeds ±1 (with 
multiphcity two), ±^j2/a, and ±{g/^Y^^ /a. The eigen- 
vectors are complete unless the eigenvalues (3-^ ± y/2a 
and ± (5/7)^/^ coincide. That is, the scalar block 
meets the criteria for strong hyperbolicity as long as 
2a^{ghY'\ 

The GH system with moving puncture gauge condi- 
tions is strongly hyperbolic everywhere, except for re- 
gions of spacetime in which 2a = {g/^y^^ ■ This restric- 
tion on strong hyperbolicity also applies to BSSN with 
the moving puncture gauge (32j . In fact, the characteris- 
tic speeds for GH with moving puncture gauge are pre- 
cisely the same as for BSSN with moving puncture gauge. 
It is recognized from studies with the BSSN equations 
that the condition 2a ^ [g/lY^^ is typically violated 
in black hole spacetimes on a 2-dimensional surface in 
space [m, [l^l . The breakdown of strong hyperbolicity 
does not appear to cause problems for finite difference 
codes. On the other hand, the lack of hyperbolicity can 
create serious problems for spectral codes that rely on 
the passing of characteristic information between spatial 
domains [33|. 

Recall that the moving puncture gauge conditions pip 
can be modified by using the constraint = to replace 

with AT'jkg^''- With this replacement Eq. jBlJl) be- 
comes 



dtP^ = -(.9/7)^/=^ 



g'^g'^'dkge, - ^g^'g^'d^gu 



(B6) 



The principal symbol (jB2p along with its scalar and vec- 
tor blocks are modified accordingly. However, the eigen- 
values arc not changed, and once again the eigenvectors 
are complete if 2a ^ {g/lY^'^- 
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